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Abstract 

We present a measurement of the branching fraction for the doubly-Cabibbo-suppressed decay 
D + — > K + ir°, using 281 pb _1 of data accumulated with the CLEO-c detector on the ^>(3770) 
resonance. We find B (D + -> K + tt°) = (2.25±0.36±0.15±0.07) x 10~ 4 , where the first uncertainty 
is statistical, the second is systematic, and the last error is due to the uncertainty in the reference 
mode. The results presented in this document are preliminary. 



* Submitted to the 33 rd International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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The Cab ibbo- favored hadronic decays of the c quark proceed through c — > slVy , Wy — > 
ltd. The doubly- suppressed decays proceed through c — > dWy, Wy — > us, and are expected 
to be suppressed by a factor \(V c dV us ) / {V cs V u d)\ 2 ~ 2.5 x 1CT 3 . We have measured the 
branching fraction for the doubly- Cabibbo-suppressed decay D + — > fC + 7r° (charge conjugate 
mode .D - — > K~tt° implied also, except where noted). 

For this measurement we have used a sample of 281 pb _1 e + e~ — > -0(3770) events, pro- 
duced with the CESR-c storage ring and detected with the CLEO-c detector. The data sam- 
ple contains about 0.8 x 10 6 D + D~ events (our target sample), one million D°D events, 
five million e + e~ — ► uu, dd, or ss continuum events, one million e + e~ — > r + r _ events, 
one million e + e~ — > jip' radiative return events (sources of background), as well as Bhabha 
events, //-pair events, 77 events (useful for luminosity determination and resolution studies). 

The CLEO-c detector is a general purpose solenoidal detector which includes a tracking 
system for measuring momenta and specific ionization (dE/dx) of charged particles, a Ring 
Imaging Cherenkov detector (RICH) to aid in particle identification, and a Csl calorime- 
ter for detection of electromagnetic showers. The CLEO-c detector is described in detail 
elsewhere jl 3HII. _ 

The resonance ■0(3770) is below the threshold for DDn production, and so the events of 
interest, e + e~ — > 0(3770) — > DD, have D mesons with energy equal to the beam energy, 
and a unique momentum. Having picked the particles being considered to make up a D 
meson, following Mark III we define the two variables: 



AE = ^Ei — -Ebeam , (1) 

i 

and 

M bc = ^ b 2 cam -|E^I 2 . (2) 

where E{, Pi are the energy and momentum of each particle making up the D meson. For a 
correct combination of particles, AE will be consistent with zero, and the beam-constrained 
mass Mb c will be consistent with the D mass. 

In addition to D + — > K + n Q , we have studied the singly- Cabibbo-suppressed decay 
D + — ► 7r + 7T°, as a higher rate decay possessing many of the features of D + — > fT + 7r°, 
and the Cabibbo-favored decay D + —* K^-k + -k + , as a high rate, clean mode used for nor- 
malization. We distinguish between K ± and tt^ using information from the RICH, and 
dE/dx information from the central drift chamber. We detect 7r°'s via tt° — > 77, detecting 
the 7 rays in the Csl calorimeter. We require that the calorimeter clusters be above 30 MeV, 
have a lateral distribution consistent with that from 7 rays, and not be matched to charged 
tracks. We require that the 77 invariant mass be within 3 standard deviations of the ir° 
mass. 

We select candidate combinations that have AE between —40 MeV and +35 MeV, for 
K + 7T° and 7r + 7r°, and between —20 MeV and +20 MeV for K~tt + tt + . These are roughly 
3 standard deviations requirements, based on Monte Carlo studies. The asymmetric cut 
for K + tt° and 7r + 7r° is due to a low-side tail on n° energies, and the wider window due to 
poorer energy resolution. To study background, we select combinations with AE between 
-100 and -50 MeV, and between +45 and +100 MeV (+50 and +100 MeV for K-n + n + ). 
On the rare occasion when an event contains more than one K + tc° combination that passes 
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our AE requirement, we choose the combination with AE value closest to zero. Multiple 
candidates per event for 7r + 7r° and for K~7r + n + are removed by the same procedure. Thus, 
we allow only one candidate per event per decay mode per D charge. 



|M BC Distribulion] M. : Distribution M 5; Distribution 
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FIG. 1: M bc distributions of D+ -> K-ir + ir+, D+ -> vr+vr and L>+ if+vr . The points are 
obtained by selecting the AE signal region, the shaded histogram is from the AE sidebands, and 
the line is the fit described in the text. 

The Mb c distributions for candidate combinations are shown in Fig. ^ The normalization 
mode D + — > K~h + 'k + is essentially background free. The D + — > 7t + tt° mode background 
is well described by the distribution obtained from the AE sideband, as is that for the 
D + — > K + n° mode. There is a clear peak in D + — > K + tt°. 

We perform an unbinned maximum likelihood fit to extract signal yields from the Mb c 
distributions. For the background, we use an ARGUS function p, with shape parameter 
determined from the AE sideband M^ c plot. For the signal, we use a Crystal Ball line 
shape 0, which is a Gaussian with a high-side tail. As Monte Carlo studies show that 
D + — > i^ + 7r° and D + — > 7i + tt° have the same signal shapes, we have determined the line 
shape parameters (Gaussian peak location, Gaussian width, point at which high-side tail 
begins) from the D + — > 7r + 7r° Mbc distribution, and used them in the fit to the D + — > K + tt° 
Md C distribution. We have varied the shape of the high-side tail, as part of the systematic 
error study. 

TABLE I: The MC efficiencies, fit yields and branching fractions. A correction has been applied 
to the branching fraction for vr°-fmding efficiency (see below). Only statistical uncertainties are 
included. 



Mode 




Yield 


B (%) 


D+ -» K-1T+7T+ 


54.91 ± 0.17 


79612 ± 290 


Input 


D+ -> vr+vr° 


50.16 ± 0.16 


963 ± 54 


0.1311 ± 0.0074 




44.53 ± 0.15 


147 ± 23 


0.0225 ± 0.0036 



Results of the fits are shown in Table HI Also given in Table |J is the detection efficiency 
for each mode, and the branching fractions obtained for D + — > 7r + 7r° and D + — > K + n Q . 
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Those branching fractions are obtained relative to D + — > K~tt + tt + , taking that branching 
fraction as (9.40 ± 0.30)%, which is obtained from a weighted average of the Particle Data 
Group (PDG) value 11 1 and the recent CLEO measurement j8[. The branching fraction for 
D + — > 7T + 7r° is in good agreement with our previously-published branching fraction using 
the same data set ((0.125±0.006±0.007±0.004)%) @, is not independent of it, and should 
not be used in place of it. 

We have considered many sources of systematic error to the D + — > K + tt° branching frac- 
tion, including: signal Monte Carlo statistics, track-finding efficiency, 7r°-finding efficiency, 
particle identification, the AE requirement, final state radiation, and the uncertainty from 
our fitting procedure (background shape, signal shape). The only ones greater than 1/10 
the statistical error are 7r°-finding efficiency, background shape, and signal shape. 

The Monte Carlo simulation of the calorimeter response to photons is imperfect, partic- 
ularly in those angular regions where there is considerable material between the interaction 
point and the calorimeter. Consequently, the Monte Carlo simulation slightly overestimates 
the efficiency for detecting 7r°'s. Various data-Monte Carlo comparisons suggest a correction 
factor of (0.95 ± 0.04), which we apply. 

The background shape is determined by a fit to the AE sideband data. The error on the 
shape parameter thus determined translates into a ±4.4% relative error in the D + — > K + tt° 
branching fraction. The signal shape is determined by a fit to the D + — > 7i + tt° signal. 
Uncertainty comes from the determination of Gaussian width (a), and point at which non- 
Gaussian tail sets in (a). We have determined the error ellipse in the determination of these 
two parameters, and noted the variation in fitted D + — > K + tt° yield as one travels around 
this error ellipse. In that way, we obtain a systematic error of ±2.6%, relative. Note that 
both the background shape uncertainty and signal shape uncertainty are really statistical 
errors, hence can be reliably determined and will decrease as additional data are taken. 

Our final result is 



B (D + -> Ji+vr ) = (2.25 ± 0.36 ± 0.15 ± 0.07) x 10~ 4 

where the first error is statistical, the second error is systematic, and the third error is from 
the uncertainty in the D + — > K~7t + tt + branching fraction, (9.40 ± 0.30)%, used as the 
normalizing mode. 

Our result is in good agreement with the only other measurement of this branching 
fraction, Babar's recent B (D + K + n°) = (2.52 ± 0.47 ± 0.25 ± 0.08) x 10~ 4 [10]. It can 
be converted to a width, using the PDG value for the D + lifetime, and compared with 
the width for doubly- Cabibbo-suppressed D° decay D° — > K + tt~ , using the PDG value for 
D° — > K + 7i~ branching fraction and D° lifetime [11]. In this way we obtain 



T(D° -> K+7T-) B(D° -> K+7i-) x r D+ 

Implications of our result are discussed in [l2j. 
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